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The use of qubits as sensitive magnetometers has been studied theoretically and recent demon- 
strated experimentally. In this paper we propose a generalisation of this concept, where a scanning 
two-state quantum system is used to probe the subtle effects of decoherence (as well as its surround- 
QQ , ing electromagnetic environment). Mapping both the Hamiltonian and decoherence properties of a 

• qubit simultaneously, provides a unique image of the magnetic (or electric) field properties at the 

' nanoscale. The resulting images are sensitive to the temporal as well as spatial variation in the fields 

^N) , created by the sample. As an example we theoretically study two applications of this technology; 

►-^ . one from condensed matter physics, the other biophysics. The individual components required to 

Q ' realise the simplest version of this device (characterisation and measurement of qubits, nanoscale 

, positioning) have already been demonstrated experimentally. 

^ ■ PACS numbers: 68.37.-d,03.65.Wj,03.65.Yz,05.40.-a 

^ ! 

,— The ability to image the very small is a powerful tool in the physical and life sciences. For each advance in imaging 
' techniques, there have been equally important steps in understanding the systems they image. We theoretically show 
how a scanning two-state quantum system provides a unique perspective through monitoring of the subtle effects of 
> ■ decoherence. Previous work has focused on using a coherent two-state probe as a sensitive spatial electrometer or 
^ ' magnetometer. The system proposed here greatly expands this view by simultaneously imaging spatial information 
and the fluctuations occurring within. To illustrate the technique, we investigate two contrasting applications - one 
^H' from the realm of condensed matter physics, the other biological. Our results demonstrate that the technique provides 
a fundamentally new nano-imaging mode with wide applicability, combining both structure and process information 
in a way not accessible by other means and, given the recent demonstration of scanning qubit magnetometry, that 
implementation is well within reach of current technology. 
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I. DECOHERENCE MICROSCOPY 



Quantum coherence and entanglement in mesoscopic systems are such delicate effects that we usually treat the 
environment as an inevitable influence to be mitigated against at all costs. One seeks to isolate the quantum system 

00 , as far as possible from the environment in order to maintain and observe non-classical effects. The imaging technique 
we introduce here is based on the inversion of this viewpoint - the inherently quantum properties of a probe system 
are deliberately exposed to an environment in a controlled fashion. The effects of the environment are monitored and 

1 J ■ analyzed as a function of position. The resulting scan forms a unique image of the position and behavior of fluctuations 
rS I in the environment as measured by the direct electromagnetic effects on the probe and its loss of quantum coherence. 

■ The use of qubits as sensitive magnetometers (or electrometers) has previously been discussedii^ii^i^i^ as the 
energy levels, and hence coherent evolution, depends very precisely on its electromagnetic environment. Of particular 
interest to our proposal is the recent demonstration of a room temperature diamond colour-centre qubit scanning 
magnetometer systemii^. Here we focus on a novel imaging mode, rather than increased sensitivity. Monitoring the 
position dependent probe Hamiltonian and decoherence channels allows one to use the scanning quantum system as 
both a magnetometer/electrometer and as a probe of the decoherence environment simultaneously. As we will show, 
this imaging mode can reveal new features not observable by other means, as it is sensitive to the temporal as well 
as spatial variation in the fields created by the sample. For a given sample, this provides more than just structural 
information, it also provides unique information on processes occurring in the sample and where these processes occur. 

The paper is organised as follows: we first present the concept of scanning quantum decoherence microscopy in a 
non-specific way and illustrate its operation with a general example. Following this, we present an example using 
existing technology, where we study a qubit composed of a nitrogen-vacancy defect within a diamond substrate. We 
simulate the use of this system in imaging a biological macromoleculc with non-trivial magnetic properties. 

The imaging setup we envisage employs a two-state system (a qubit), based on charge or spin degrees of freedom, 
attached to the end of a cantilever probe, as illustrated in Fig.[Ua). The quantum probe is governed by the Hamiltonian 
Hp = ea^+^cTz, where the a's are the Pauli operators and the coefficients (e and A) are functions of both the operating 
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FIG. 1: Schematic of a scanning decoherence microscopy setup, (a) The probe qubit is scanned across the sample while its 
quantum state is (weakly) monitored, (b) At each point {x,y) a measurement record Ixy{t) is obtained. Using the time- 
correlated signal, the spectral response Sxy{io) of the qubit probe is determined, (c) From this data, a measurement of the 
effective qubit Hamiltonian Hxy (static magnetometer or electrometer) field map and decoherence rate V^y map as a function 
of probe position is obtained. In this figure, the sample consists of fiuctuators with a non uniform distribution in both spatial 
density and fluctuation rate. The decoherence map reveals the fluctuator distribution not visible in the static field map. (d) 
Combining this information provides a direct window into the distribution and character of the sources of field fluctuations in 
both space and time. In this image, the vertical scale is given by the strength of the field, whereas the colouration is given by 
the effective decoherence rate. 



conditions and the probe's electromagnetic environment. The probe is scanned across a sample and at each point, the 
Hamiltonian and decoherence felt by the qubit is determined through a series of measurements. The decoherence of 
the qubit, in particular, is highly sensitive to the fluctuations in the local environment. Thanks (in part) to the ongoing 
work to construct controllable quantum devices^, a number of techniques have been developed to measure the state or 
evolution characteristics of an (open) quantum system. These include techniques to reconstruct an arbitrary quantum 
state or process^iiiS, quantum channel^iii and the Hamiltonian governing a few state quantum syste m^^i-'^^i-'^^i-'^^i-'^^i-'^''' 
or even spin echo techniques'^ from magnetic resonance. 

The simplest realisation of decoherence microscopy is to perform a Rabi oscillation or spin echo experiment on a qubit 
which can be moved relative to the background decoherence source. In this case, the measurement procedure (initialise, 
evolve, measure) is used to map the (ensemble averaged) time evolution of the system and the decoherence rate is fitted 
as an exponential decay. With sophisticated measurement schemes and fitting procedures, the image can be acquired 
more quickly and efficiently. In this paper, we will focus on the technique of Hamiltonian characterisationi2ii2iiil and a 
general weak continuous measurement of the qubit as an example. It is important to note that the particular choice of 
continuous measurement and Hamiltonian characterisation is not a requirement to perform decoherence microscopy, 
it simply provides an elegant and efficient realisation but other more established experimental techniques could be 
similarly applied. 

We consider a quantum probe, located at position {x,y), suspended a distance hp above the sample. The probe 
is weakly monitored, providing a measurement record 2xy(t). This measurement record is then correlated and the 
spectral response Sxy{i^) computed, see Fig.[T]Jb). The Hamiltonian and decoherence parameters can then be estimated 
from this spectral response, for each probe position. If the components of the Hamiltonian are plotted as a function 
of position across the sample (Hxy), the probe acts as a sensitive electrometer or magnetometer (depending on the 
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type of qubit). as given in the H^y example in Fig. [TJc). However, complete analysis of the probe evolution allows the 
decoherence channel(s) to be extracted (T^y), giving information about the strength, direction and character of the 
dynamics of the environment, as well as the induced static field. The resulting decoherence scan, shown in Fig.[ljc), 
reveals new information about the fluctuator frequency distribution in the sample which was not apparent in the 
electrometer/magnetometer image. Combining both sets of data, Fig. [ijd), provides a direct image of the correlation 
between the temporal dynamics of the environment and the spatial structure. 

To make our proposal quantitative, we introduce a model for weak continuous measuremeni^iS which captures all 
of the essential physics*. Here, we assume the measurement of the qubit can be modeled as a inefficient (or weak) 
positive operator valued measure (POVM) in the (7z basis. The density matrix after measurement, p' , is given by 



Tr[AlA±p] 



t 

(1) 



for a measurement operator 



A± = ^(^/^±7^|0)(0| + ^/T^|l)(l|), (2) 



with some measurement strength k. The measurement process consists of repeated weak POVM measurements 
separated by a time interval At, during which time the system undergoes normal evolution. The measurement 
repetition interval is then a measure of the bandwidth of the detector, BW = 1/At. 

The measurement record I{t), is the result, +1 or —1, of a measurement at time t. The steady-state autocorrelation 
of this measurement signal X{t) is then given by 

(X(t)X(i + t)),, = Tr[CT,e^"(7zP.s] (3) 

via the quantum regression theorem^Si^. Here pss = p{oo) is the steady state density matrix and e^"^ is the solution 
to the density matrix evolution governed by p = C[p,H,T], where H is the qubit Hamiltonian and F represents the 
decoherence rate of the system without measurement. The spectrum of the signal is then 



S(uj) J' 



n{<yAt))l (4) 



where J-'[{az{t))] is the Fourier transform of the (ensemble averaged) expectation value of the ct^ operator. From this 
response spectrum, we extract the Hamiltonian and decoherence parameters directlyi^ii^ for each spatial location 
across the sample. When performing this parameter estimation process, it is not necessary to make any assumptions 
about the source or mechanisms producing the decoherence. Simply mapping the induced decoherence as a function 
of position provides a dynamical image of the sample fluctuations. 

In the limit of small k, the model of Eq. ([T]) is equivalents^ to more complicated master equation mode b^^i^^'^'^i^^'^^ . 
Expanding the evolution to first order in both k and At, we can derive an equivalent Lindbladian master equation 
with an effective CTz decoherence channel of strength 

^2 

Fmcas ~ 4^^' ^^-^ 

which corresponds to the measurement induced decoherence. The measurement strength should be chosen such that 
this is smaller than both the sample induced decoherence and other intrinsic decoherence sources^. The intrinsic noise 
sources include those introduced by the tip supporting the qubit, although ideally this decoherence is equal to or less 
than that introduced by a bulk substrate, owing to the geometry of the tip. 

We can also calculate the information extracted at each measurement step. At, by looking at the reduction in 
entropy, AS'e, of the system. Expanding for small k, this gives 

A5e = 5e(p) - 5e(p') = + 0{k^) (6) 

loge(4) 



* While this is the more general scenario, for some systems strong projective measurement is a more natural implementation. In this case, 

projective or 'DC measurements can be used to measure the system evolutioni^ii^iiS, producing equivalent results, 
t This means operating well above the Korotkov-Averin boundSSiS so that the sample decoherence effects dominate the signal. 
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as the information obtained (in bits) from a single measurement of an initially mixed state p = (|0)(0| + |l)(l|)/2. As 
we increase k or the bandwidth, the amount of information obtained in a given time interval is increased, Eq. ([S]), at 
the expense of greater measurement induced decoherence, Eq. ([5]). 

The spatial resolution of the probe system is ultimately governed by the effective strength of the environmental 
decoherence as a function of distance. Most decoherence channels (for solid-state qubits) depend, in some way, on the 
inverse of the separation between qubit and environment. The coupling between a single decoherence (point) source 
locate at Xs and the probe qubit is, in general, proportional to l/r", hence the response of the qubit as a function of 
(a;-)position across the sample is given by a exponentiated Lorenztian, 



nix) 



1 



+ (x - XsY 



i/2 



(7) 



The Full Width at Half Maximum (FWHM) of this function gives the spatial resolution Ax = 2/ip\/2^/" — 1, where 
hp is the height of the probe above the sample. By inspection, we see an electric dipole induced potential [l/r^) has 
a FWHM Ax = 2hp whereas a magnetic dipolar interaction has slightly better resolution. Ax l.bihp. 

We now consider two examples in which we simulate the effects of a decohering environment on an example probe. 
These examples provide both a straightforward illustration of the power of imaging the sample induced decoherence 
and a test of the feasibility using current and near future technology. 



II. EXAMPLE I: IMAGING THE DISTRIBUTION OF BACKGROUND CHARGE FLUCTUATORS 

Our first example system comprises an electrostatic qubit interacting with a sample containing background charge 
fluctuators. This example is of particular interest as background charge fluctuations have been the subject of intense 
scrutiny due to their relevance to solid-state quantum devices^^i^S. In this case, we can probe the spatial and frequency 
distribution of these fluctuators in a way which is not possible using current microscopy techniques. 

A suitable probe qubit would be any of the myriad of charge based qubit designs such as quantum-dots^^, donors^ 
or coopcr-pair-boxs^i^i^. A probe based on a cooper-pair-box (CPB) system provides a particularly good example 
as CPB qubits are now regularly produced exp er iment all y ' ^ ^ i '^^ and the bias point of the system can be varied, 
resulting in a change in the sensitivity to different components of the environmental decoherence^l. For generality, 
we model the qubit purely as a two-state system which interacts via a field dependent component in its Hamiltonian 
(see appendix]^ and therefore all the dimensions in this example are given in terms of a normalized length scale (L). 

Let us consider a 2D surface which contains background charge fluctuators that we wish to image. We use a 
simplified model of these fluctuators (see appendix |^ in the fast fluctuator limit, to illustrate the concepts. In 
example II, we will consider a more specific example in the slow fluctuator limit. We take the potential felt by the 
qubit as a simply electrostatic potential due to a charge defect dipole and simulate the effective decoherence field felt 
by the qubit due to the fluctuator bath (defects in the sample). In Fig. [2{a), we have an fictitious sample comprised 
of regions (outlined) containing a uniform fluctuator distribution with an area density of 1000 defects per square (all 
distances are in normalised units). The frequencies of these fluctuators are then assumed to be distributed with a 1// 
distribution. We have also included a 'calibration' region on the left of the sample, where the frequency of fluctuators 
is varied in a controlled fashion. This will illustrate the frequency selective nature of the measurement process. 

The variation of the probe Hamiltonian, H^y^ due to the presence of background charges is shown in Fig. [U^b). 
The background charges are taken to be charge dipoles^Sii^i coupling to A, the az component of the qubit, as l/r^ 
(see appendix El and the probe is positioned hp = 0.05L above the sample. Initially assuming the fluctuators to be 
static, the spatial variation in the Hamiltonian coupling term is then a measure of the electric potential induced by 
the dipolcs (a simple electrometer). The state of the fluctuators will not be static (in general) and the total induced 
field will result from an average over the fluctuator states. In Fig.[2Kc), we calculate the decoherence felt by the probe 
due to the combined effects of all the fluctuators in our fictitious sample. The resolution for this example increases by 
a factor of 1.55, compared to measuring the Hamiltonian, as F cx S{uj) oc 1/r^. Comparing these images, we see that 
the decoherence measurement is more sensitive to the fluctuators whose frequency arc closest to the probe energy. It 
is known that, for a 1// bath, a relatively small fraction of the total fluctuators contributed a large amount of the 
total decoherence^l but here we see it strikingly depicted in the images. 

To further demonstrate the utility of mapping both Hamiltonian and decoherence components simultaneously. 
Fig. [2jd) is a combined plot of both Fig. [2jb) and Fig. ^c). In this case, the effective field (static) is a surface plot 
with colouration corresponding to the decoherence rate (dynamics). The vertical scale of this 'mountain range' plot 
corresponds to the varying electric potential across the sample. The red 'mountains' are fluctuators with frequency 
close to the qubit transition frequency while blue corresponds to far off resonant fluctuators. In principle, the frequency 
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FIG. 2: (a) A fictitious sample with uniform distribution of fluctuators (coloured dots) within non-trivial spatial regions 
(outlined in white). The fluctuators are modeled as point sources but have been enlarged here for clarity. As we have not 
chosen a particular qubit architecture, all distances are in normalized length units (L) and the colour coding of the fluctuator 
locations depends on the frequency of the fluctuator, in units of the qubit probe frequency. The region on the left is a 
calibration region where the fluctuator frequency varies with position in the region. The other regions contain fluctuators with 
a uniform spatial and 1// frequency distribution, (b) Simulated response of the probe Hamiltonian interacting with a uniformly 
distributed bath of charge dipoles, with hp = 0.05L. The normalised variation of the Hamiltonian component (A) induced by 
the background charges is plotted as the colour scale, (c) The effective decoherence (normalized) rate felt by the qubit as a 
function of position. The qubit is now sensitive to the 1// nature of the fluctuator bath, with a relatively small number of 
fluctuators contributed a large fraction of the decohering effect, (d) Plots (b) and (c) have been combined using the value of 
(b) for the height and (c) for the colour. The frequency selective nature of the image is apparent, where fluctuators closest 
to the probe in energy contribute the most decoherence signal (red peaks) where other peaks (blue) are to due to fluctuators 
which do not strongly decohere the qubit. 

sensitivity of the qubit can be tuned to probe different components of the decoherence, to obtain more information 
about the spectral characteristics (1//, Ohmic etc. ) as a function of position. 

III. EXAMPLE II: IMAGING THE POSITION AND SPIN OCCUPATION OF FERRITIN MOLECULES 

Now, we choose to use both a probe qubit and sample which have been well studied experimentally. In contrast 
to the previous example, where all parameters were dimensionless, here we investigate a specific implementation to 
evaluate the technical feasibility of the proposal. We consider the imaging of (bio)molecules with large uncompensated 
spin, such as Horse-spleen Ferriti n^^'^^ or Feg^. Here, the point is neither to image an individual spin^ or image the 
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FIG. 3: Variation in strength of the maximum signal {<j} = 0) as a function of probe height for different sized spins. The 
signal is the fractional variation in the probe spin Hamiltonian due to the influence of the sample spin. The upper and lower 
threshold limits come from the measurement bandwidth and decoherence respectively. The insert illustrates the simplified 
magnetic dipolar coupling model, which depends on the separation, orientation and magnetic moment of both the probe and 
sample spins. 



location of the molccule a^^i'^^' as both can be done with existing technology. We show that qubit probe imaging can 
both map the location of the spins and probe their magnetic dynamics. 

The decoherence introduced at the probe qubit will be a function of both the interaction strength and the flipping 
rate of the sample spin(s). As wc are considering large sample spins in a static magnetic field at low temperature, 
we will assume that the flipping rate is slow on the scale of the probe Hamiltonian. The spectral response of the 
qubit is split, with the separation between the peaks giving the effective difference in the Hamiltonian between the 
two sample spin states (see appendix [B]) . We compute the coupling strength between the spins, given the simplified 
magnetic dipolar model, as illustrated in the insert of Fig. [3) The ability to resolve the induced coupling is ultimately 
limited by the total effective decoherence rate of the qubit. The ratio of the peaks also gives a measure of the relative 
spin populations, which in turn relates to the effective temperature and/or flipping mechanism. 

We use a nitrogen-vacancy (NV) centre in diamond as our probe spin, as this has been shown to be a controllable, 
well isolate spin system which displays stable quantum coherent properties up to room temperatur o^^i^^i^^i'^'''i^^i^^i^°i^-'^ 
and has even been used as a nanoscale scanning magnetometer^^, in a similar cantilever configuration to that which 
we envisage. We then couple this probe to a bulk spin of order Mq = 50 — 200yUB and include the effects of intrinsic 
decoherence and finite measurement bandwidth. Using experimentally realistic parameters for both the probe qubit 
and the sample spin, we plot the response of the system as it passes over the spin. An alternative probe qubit system 
would be to use a microSQUID or flux qubili^i^i^. 

For our example, we use known system parameters for an NV centre driven by a microwave loop and readout via a 
laser probe measurement^. Wc take the Rabi frequency of the qubit to be lOMHz and the measurement bandwidth 
BW = lOOMHz but not necessarily 'strong'. The intrinsic decoherence rate is approximately 100 times slower than 
the Rabi frequency. The measurement strength k is chosen such that the measurement induced decoherence is weaker 
than the intrinsic decoherence, for a given detector bandwidth. 

Fig. [3] shows the maximum fractional Hamiltonian component as a function of probe height for three different spin 
samples. The measurement window is defined between the decoherence rate and the measurement bandwidth. The 
probe height can vary over almost 100 nm and still provide a detectable signal. In Fig. Ul we create an image from a 
fictitious sample consisting of 4 mesoscopic spins with varying net magnetisation. Wc assume the spins are in thermal 
equilibrium with the sample substrate (which we have set at T = AK) and that they arc fiipping due to thermal 
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FIG. 4: (a) Magnetic field strength as measured via the variation in the probe qubit Hamiltonian. The pixel intensity represents 
the normalized fractional variation in the probe spin Hamiltonian. The detection limits discussed earlier for maximum and 
minimum detectable field are included in this calculation, (b) The effective temperature of each spin, based on the population 
of the ground and excited states. The intensity in this plot is given by the effective field strength of each spin signal separately. 
The white/black circles indicate the position and diameter (8nm) of the sample spin and enclosing molecules, (c) Surface plot 
of the magnetisation with colouration based on the fraction of time spent in the excited spin state. The height and colourscale 
are normalised as per (a) and (b) respectively. The grey colouration corresponds to a region where the peak splitting is not 
large enough to resolve the spin populations with these measurement parameters. 

processes-'-. The magnetisations arc Mq = 50, 70, 100 and 200/iB and the average population of the excited state is 
given by a Bohzmann distribution for a background magnetic field of B = O.IT and temperature T = 4K. The spatial 
resolution of the probe position is a 50 x 50 grid, giving 2500 points over lOOOOnm^ and the probe height -was set to 
hp = 20nm. Fig. Ufa) sho-ws the measured magnetic field over the sample. Note that the probe in this mode (purely 
acting as a magnetometer) does not successfully resolve two of the spins. 

Fig. sho-ws the measured decohcrence field over the same sample. As each spin has a different magnetisation, 
the deeoherence effects (in this case splitting of the spectral response peak) resulting from each spin are different. The 



The exact details of the flipping process are unimportant for this example. We assume they are dominated by thermal processes and 
obey Boltzmann statistics. 
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FIG. 5: The image from Fig. HJa) with simulated noise stemming from a finite measurement time. The images from left to 
right simulate dwell times of fdwoii = 2/is, 200/is and 20ms which give a total acquisition time for a 50 x 50 grid of ta ~ 5ms, 
0.5s and 50s respectively. The average pixel intensity variance is 10"^, 10"^ and lO"'' for the three images. 



ratio of the two split peaks provides the population of the spin states, which is directly related to the magnetisation 
and effective temperature of the sample spin. In this plot, the ratio of the split peaks has been used to colour code 
the data, with blue indicating a large spin magnetisation (or low effective temperature) and red a small magnetisation 
(high temperature). The intensity of the colour is given purely by the amount of signal available from each decohcrence 
source (compared to the probe spin's intrinsic decohcrence), whereas in Fig. Sl^a), the intensity was proportional to 
the total induced field. Finally, we can combine this data to produce a plot showing the field intensity with each 
decohcrence source (mesoscopic spin) tagged based on its effective temperature. This is shown in Fig.SJc), where the 
existence of all four spins can be detected based on the colour tagging, in contrast to the magnetometer scan alone. 



IV. IMAGE ACQUISITION RATE AND NOISE 



While we have demonstrated that new information can be obtained by looking at the induced decohcrence, this 
is only useful if the information can be obtained within an experimentally accessible time. Using the measurement 
model discussed earlier, we can estimate the parameter uncertainties in the Hamiltonian characterisation proces a^'^'^^ . 
Retaining the parameters from Example II, we calculate the noise expected for a finite dwell time (tdweii) on each 
pixel and the total image acquisition time {ta). Fig. [5] shows Fig.[3]Ja) with the noise resulting from a finite bandwidth 
and dwell time. The image acquisition times (for 2500 pixels) varies between 5ms and 50s for pixel error variances 
which range from 10"^ to 10"^. 

To allow fast image acquisition times (using continuous a measurement mode), it is important to have both large 
detector bandwidth and a large ratio of bandwidth to intrinsic decohcrence rate BW/Fq. The qubit transition 
frequency is less important, provided it is at least an order of magnitude greater than the intrinsic decohcrence. 
However, a tunable qubit frequency is advantageous as many types of decohcrence depend strongly on the frequency 
at which they are probed. 

To estimate detector noise as a measure (rather than acquisition time), we take recent noise limit estimates for 
magnctometry based on NV centres^ as a guide. In this case, the measurement is a 'DC measurement (BW = lOkHz) 
with a sensitivity of 20nT/-\/Ha^ and a lOMHz qubit frequency. This results in an additional contribution to the pixel 
variance of « lO"^ but still results in a clear image (see Fig. [5]). 



V. CONCLUSIONS 



We have proposed and theoretically investigated a fundamentally new and complementary imaging mode which 
takes advantage of current qubit technology and modern techniques for characterising few-state quantum systems. 
Mapping the induced decohcrence across a sample indirectly images the dynamics of the environment, providing 
a new window into the microscopic world with wide applications to spin and charge fluctuations and transport in 
both solid-state and biological systems. In contrast to a scanning qubit magnetometer, the decoherence signature is 
sensitive to fluctuations even when the time-averaged field is close to zero. Such processes abound in nature, in which 
the time dependent behaviour of electromagnetic fields stems directly from dynamical processes of interest, such as 
ion flow, spin fluctuations or conformational change. The device may be realised using a range of architectures, each 
of which will be sensitive to a particular electromagnetic effect within a sample. While we have discussed several 
spcciflc examples using near future technology, the basic requirements for a decoherence microscopy already exist 
in the form of nano-scale positioning and individual two-state quantum systems which can be initialized, controlled 
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and measured. Hence, this technique has the potential to open a new window for imaging nanoscale processes in the 
physical and life sciences. 

Acknowledgments 

The authors would like to thank N. Oxtoby, C. Hill, G. Milburn and G. Schon for useful comments on the manuscript. 
This work was supported by the Australian Research Council. JHC acknowledges the support of the Alexander von 
Humboldt Foundation and LCLH is the recipient of an Australian Research Council Australian Professorial Fellowship 
(DP0770715). 

APPENDIX A: EXAMPLE I - DECOHERENCE DUE TO A BATH OF 1// FLUCTUATORS 

We wish to use a simple model of the decoherence felt by a probe (charge- )qubit interacting with a bath of 1// 
fluctuators22i^i22. The numerical parameters used for such a model vary greatly depending on the system and even 
from sample to sample. Here we are interesting in the functional dependence, rather than the explicit values. We 
draw heavily from references [37j and (28| as an example. 

Consider a qubit interacting with a bath of fluctuators via an interaction Hamiltonian Hmt which has the form 

i^i„t = (T.^w,6]foj (Ai) 

where az acts on the qubit and bj destroys (creates) an electron in a localized state in the bath. 
The spectral response from the jth fluctuator (in the fast fluctuator limit) is given by 

where vj is the strength of the fluctuator which fluctuates with rate 7^. The total spectral response is then a sum 
over each of the fluctuators, 

5(^)=^s,(c.) (A3) 

j 

and a simplified Golden rule model gives the relaxation rate (F';-'^) and dephasing rate (F^''^) as 

= ir?« = \s{Ej) (A4) 

where Ej is the tunnelling energy of the qubit. Given a functional form for how Vj{r) depends on the fluctuator/qubit 
separation r, we can then compute the effective decoherence rate felt by the qubit. 

APPENDIX B: EXAMPLE II - MAGNETIC DIPOLAR COUPLING BETWEEN A PROBE AND 

SAMPLE SPIN 

For this analysis, we will use a simplified (but quite general model) consisting of a probe spin interacting via the 
magnetic dipolar interaction with a much larger sample spin (see insert of Fig. [31). The sample spin can be considered 
as a Ferritin, Fcg or other mesoscopic molecule with a net magnetic moment Mq ^ IfiB- As the spin is large (and to 
simplify the analysis of decoherence) we will assume that it is in thermal equilibrium with the sample environment 
and therefore the average magnetisation and spin flip rates are given by the standard thermodynamic quantities. In 
making this assumption, we ignore quantum mechanical effects between the probe and sample spins and treat the 
effect of the sample spin on the probe as a fluctuating classical field. 

The magnetic dipolar interaction is given hy^ 



P.S 3{P.f){S.r) 



(Bl) 
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where 7p and 7s are the probe and sample spin gyromagnetic ratios, r is the vector separation between the spins and 
P and S are the probe and sample spin operators. 

We assume the system is bathed in a global magnetic field B which orientates both the sample and probe spins^ 
and sets their energy scales. The dipolar interaction is then given by the separation between the spins r and the angle 
subtended between the spin orientation and the vector separating the spins 4>, as illustrated in the insert to Fig. [31 
The simplified coupling Hamiltonian is therefore 

^-p=O(^)(3cos^0-l)P..5. (B2) 

which is a purely Ising type interaction whose strength depends on both the separation and angle between the spins. 
The perturbing effect of this dipolar Hamiltonian can then be measured in the spectral response of the qubit, giving 
a direct link to both the spin state and magnetisation of the sample spin. 



REFERENCES 



1 B. M. Chernobrod and G. P. Berman, J. Appl. Phys. 97, 014903 (pages 3) (2005). 

2 E. Il'ichev and Y. S. Greenberg, Euro. Phys. Lett. 77, 58005 (5pp) (2007). 

3 V. I. Shnyrkov and S. I. Melnik, Low Temp. Phys. 33, 15 (2007). 

* R. DeSousa, arXiv:cond-mat/0610716 (2006). 

5 C. L. Degen, Appl. Phys. Lett. 92, 243111 (2008). 

^ J. M. Taylor, P. Gappellaro, L. Childress, L. Jiang, D. Budker, P. R. Hemmer, A. Yacoby, R. Walsworth, and M. D. Lukin, 
Nature Physics, 4, 810 - 816, (2008). 

G. Balasubramanian, I. Y. Chan, R. Kolesov, M. Al-Hmoud, J. Tisler, C. Shin, C. Kim, A. Wojcik, P. R. Hemmer, A. 
Kriiger, T. Hanke, A. Leitenstorfer, R. Bratschitsch, F. Jelezko, and J. Wrachtrup, Nature, 455, 648-651, (2008). 

* J. R. Maze, P. L. Stanwix, J. S. Hodges, S. Hong, J. M. Taylor, P. Gappellaro, L. Jiang, M. V. Gurudev Dutt, E. Togan, 
A. S. Zibrov, A. Yacoby, R. L. Walsworth, and M. D. Lukin, Nature, 455, 644-647, (2008). 

® M. A. Nielsen and I. L. Chuang, Quantum computation and quantum information (Gambridge University Press, Cambridge, 
2000). 

M. Paris and J. Rehacek, Quantum state estimation, Lecture notes in physics, 649. (Springer, Berlin ; New York, 2004). 
" M. Mohseni, A. T. Rezakhani, and A. Aspuru-Guzik, Phys. Rev. A 77, 042320 (pages 5) (2008). 
S. G. Schirmer, A. KoUi, and D. K. L. Oi, Phys. Rev. A 69, 050306(R) (2004). 

J. H. Gole, S. G. Schirmer, A. D. Greentree, G. J. Wellard, D. K. L. Oi, and L. G. L. Hollenberg, Phys. Rev. A 71, 062312 

(2005) . 

J. H. Gole, A. D. Greentree, D. K. L. Oi, S. G. Schirmer, G. J. Wellard, and L. G. L. Hollenberg, Phys. Rev. A 73, 062333 

(2006) . 

J. H. Gole, S. J. Devitt, and L. G. L. Hollenberg, J. Phys. A 39, 14649 (2006). 

S. J. Devitt, J. H. Gole, and L. G. L. Hollenberg, Phys. Rev. A 73, 052317 (2006). 

T. F. Jordan, A. Shaji, and E. G. G. Sudarshan, Phys. Rev. A 76, 012101 (pages 6) (2007). 

G. P. Poole and H. A. Farach, The theory of magnetic resonance (Wiley-Interscience, New York, 1972). 

T. A. Brun, Am. J. Phys. 70, 719 (2002). 

G. W. Gardiner, Quantum noise, Springer series in synergetics ; v. 56. (Springer- Verlag, Berlin ; New York, 1991). 

M. O. Scully and M. S. Zubairy, Quantum optics (Gambridge University Press, Gambridge, 2006). 
2^ A. N. Korotkov and D. V. Averin, Phys. Rev. B 64, 165310 (2001). 

Q. Zhang, R. Ruskov, and A. N. Korotkov, Phys. Rev. B 72, 245322 (pages 11) (2005). 
2"* N. P. Oxtoby, H. M. Wiseman, and H.-B. Sun, Phys. Rev. B 74, 045328 (pages 11) (2006). 
2^ N. P. Oxtoby, J. Gambetta, and H. M. Wiseman, Phys. Rev. B 77, 125304 (pages 11) (2008). 

K. S. Birdi, Scanning Probe Microscopes: Applications in Science and Technology (CRG Press, Boca Raton ; London ; New 

York ; Washington D.G., 2003). 
2^ J. Schriefi, Y. Makhlin, A. Shnirman, and G. Schon, New J. Phys. 8, 1 (2006). 

^® Y. M. Galperin, B. L. Altshuler, J. BergU, and D. V. Shantsev, Phys. Rev. Lett. 96, 097009 (pages 4) (2006). 
D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 (1998). 



^ We ignore the effects of crystal field terms on both the orientation of the sample and probe spin, though this does not decrease the 
generality of the result. 



11 



L. C. L. HoUenberg, A. S. Dzurak, C. Wellard, A. R. Hamilton, D. J. Reilly, G. J. Milburn, and R. G. Clark, Phys. Rev. B 
69, 113301 (2004). 

Y. Makhlin, G. Schon, and A. Shnirman, Rev. Mod. Phys. 73, 357 (2001). 

J. Q. You, X. Hu, S. Ashhab, and F. Nori, Phys. Rev. B 75, 140515 (pages 4) (2007). 

J. Koch, T. M. Y., J. Gambetta, A. A. Houck, D. I. Schuster, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, and R. J. 

Schoelkopf, Phys. Rev. A 76, 042319 (pages 19) (2007). 

Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Nature 398, 786 (1999). 

T. Yamamoto, Y. A. Pashkin, O. Stafiev, Y. Nakamura, and J. S. Tsai, Nature 425, 941 (2003). 

J. Majer, J. M. Chow, J. M. Gambetta, J. Koch, B. R. Johnson, J. A. Schreier, L. Frunzio, D. I. Schuster, A. A. Houck, 
A. Wallraff, et al., Nature 449, 443 (2007). 

E. Paladino, L. Faoro, G. Falci, and R. Fazio, Phys. Rev. Lett. 88, 228304 (2002). 

D. D. Awschalom, J. F. Smyth, G. Grinstein, D. P. DiVincenzo, and D. Loss, Phys. Rev. Lett. 68, 3092 (1992). 

J. Tejada, X. X. Zhang, E. del Barco, J. M. Hernandez, and E. M. Chudnovsky, Phys. Rev. Lett. 79, 1754 (1997). 

E. d. Barco, J. M. Hernandez, J. Tejada, N. Biskup, R. Achey, I. Rutel, N. Dalai, and J. Brooks, Phys. Rev. B 62, 3018 
(2000). 

D. Rugar, R. Budakian, H. J. Mamin, and C. B. W., Nature 430, 329 (2004). 

H. G. Hansma and J. H. Hoh, Ann. Rev. Biophys. Biomol. Struct. 23, 115 (1994). 

C. I. Pakes, D. P. George, S. Ramelow, A. Cimmino, D. N. Jamieson, and S. Prawer, J. Mag. Mag. Mat. 272-276 (1), 
E1231 (2004). 

F. Jelezko, C. Tietz, A. Gruber, I. Popa, A. Nizovtsev, S. Kilin, and J. Wrachtrup, Single Mol. 2, 255 (2001). 
F. Jelezko, T. Gaebel, I. Popa, A. Gruber, and J. Wrachtrup, Phys. Rev. Lett. 92, 076401 (2004). 

F. Jelezko, T. Gaebel, I. Popa, M. Domhan, A. Gruber, and J. Wrachtrup, Phys. Rev. Lett. 93, 130501 (2004). 
F. T. Charnock and T. A. Kennedy, Phys. Rev. B. 64, 041201 (2001). 

T. A. Kennedy, F. T. Charnock, J. S. Cohon, J. E. Butler, R. C. Linares, and P. J. Doering, Phys. Status Sohdi B 233, 416 
(2002). 

T. Gabel, M. Domhan, I. Popa, C. Wittmann, P. Neumann, F. Jelezko, J. R. Rabeau, N. Stavrias, A. D. Greentree, S. Prawer, 
et al., Nat. Phys. 2, 408 (2006). 

R. Hanson, F. M. Mendoza, R. J. Epstein, and D. D. Awschalom, Phys. Rev. Lett. 97, 087601 (2006). 
R. Hanson, V. V. Dobrovitski, A. E. Feiguin, O. Gywat, and D. D. Awschalom, Science 320, 352 (2008). 
J. Wrachtrup and F. Jelezko, J. Phys. Cond. Matt. 18, S807 (2006). 
A. N. Korotkov, Phys. Rev. B 63, 085312 (2001). 



